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Abstract 

We present a progress report on the Cluster Processor, a special- 
purpose computer system for the Wolff simulation of the three-dimensional 
Ising model, including an analysis of simulation results obtained thus 
far. These results allow, within narrow error margins, a determination 
of the parameters describing the phase transition of the simple-cubic 
Ising model and its universality class. For an improved determination 
of the correction-to-scaling exponent, we include Monte Carlo data for 
systems with nearest-neighbor and third-neighbor interactions in the 
analysis. 
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1 Introduction 



The goal of the Cluster Processor project is to obtain an accurate description 
of the phase transition of the three-dimensional Ising model, in particular 
to determine the universal constants that apply as well to many other phase 
transitions that occur in nature. The Cluster Processor uses special-purpose 
hardware executing the Wolff version Q of the 'cluster' Monte Carlo algo- 
rithm first introduced by Swendsen and Wang [0]. 

The Cluster Processor project was planned to be executed in two stages. 
The first stage was to construct, test, debug and operate a single processor, 
called the Prototype Processor. The second stage comprised the construc- 
tion of a system of ten of such processors that could be operated in parallel, 
which would thus allow an even higher level of statistical accuracy. 

The Monte Carlo results obtained thus far with help of the Cluster Pro- 
cessor do indeed allow a numerical determination of the magnetic and tem- 
perature renormalization exponents yu and y% with a level of accuracy ex- 
ceeding other Monte Carlo studies available to us at the time that this paper 
was submitted. Nevertheless, the results are not yet quite of the accuracy 
originally expected for the Cluster Processor. Thus, we consider the Clus- 
ter Processor project as being yet unfinished. Although the construction of 
hardware parts of this special-purpose computer system, which includes a 
total of 11 processors, has already been completed some time ago, only two 
of these were available in time to contribute to the Monte Carlo analysis in 
the present article. Two more processors have become operational in the 
meantime. 

The construction part of the first stage was completed in 1995 Q, and 
subsequent simulations on the prototype processor yielded good results ||. f|] . 
The subsequent construction of the printed-circuit boards of 10 clones of the 
prototype processor was completed in about a year and a half. However, 
during the last two years, the progress in the debugging of the hardware has 
been rather slow. This problem is caused by the absence of one of the one of 
us (A.L.T.) who built the hardware. This author is now occupied by duties 
unrelated to the Cluster Processor. 

In Section 2 we present a finite-size analysis of the Binder cumulant [^] 
of the 3-D Ising model with nearest-neighbor interactions K nn only. The 
data sample consists of results taken with the Cluster Processor, mostly 
for finite sizes L = 128 and 256, supplemented by data for relatively small 
system sizes, computed on workstations and personal computers. In Sec- 
tion 3 we extend the parameter space to include third-neighbor interactions 
i^3 n which are known to reduce the corrections to scaling. The numeri- 
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cal data for these systems are restricted to L < 128 and were obtained for 
K^, n /K nn = 0.1, 0.2, 0.3 and 0.4, using workstations and personal computers. 
The Cluster Processor is restricted to simulations of the nearest-neighbor 
model with finite sizes L = 16, 32, 64, 128 and 256. These additional 
data enable us to perform an improved analysis, in particular concerning 
the accuracy of the correction-to-scaling exponent. In Section 4 we present 
analyses of the magnetic susceptibility and of the temperature-derivative of 
the Binder cumulant [jEJ, leading to new estimates of the renormalization 
exponents yh and y%. Section 5 discusses and summarizes our findings; in 
particular see Table [/]. 

The analysis presented in Section 3 is related to, and can be seen as an 
extension of the analysis given in Ref. to a much larger body of Monte 
Carlo data. In the latter paper it was demonstrated that the amplitude of 
the leading correction to scaling can be tuned by modifying the Hamilto- 
nian: in general, the use of a two-parameter Ising-like Hamiltonian enables 
one to suppress both the temperature field and the leading irrelevant field. 
In particular the latter condition can only be realized with a limited numer- 
ical accuracy, so that it remains necessary to include the amplitude of the 
leading correction to scaling in the analysis. As emphasized in Ref. ||, the 
importance of reducing the irrelevant field lies in the fact that the results of 
the fitting procedures then become almost independent of the precise value 
of the irrelevant exponent yi, which is known only with a modest precision. 

After the work of Ref. || , several new analyses of the 3-D Ising model 
and its universality class have been performed. A modified Monte Carlo 
renormalization analysis ||, which was set up such as to suppress correc- 
tions to scaling, both at the simulated Hamiltonian and at the fixed point of 
the transformation, yielded results for the critical exponents with a precision 
comparable to those of Ref. [|[. An analysis of the spontaneous magneti- 
zation and other data obtained by the prototype Cluster processor led to 
somewhat better accuracies. Very recently, two new Monte Carlo analyses 
have been reported. Ballesteros et al. j?]] obtained very accurate numerical 
data for the nearest-neighbor Ising model, and their results do indeed have 
small error margins. Hasenbusch et al. || also simulated the two-parameter 
Hamiltonians of Ref. ]^| , using a new method to determine the parameters 
such as to suppress the leading irrelevant field. They reported results with 
a precision comparable to that of Ref. 0. A further discussion of these new 
analyses, and their similarities and dissimilarities with the present analysis 
will be included in Section 4. 

Further recent developments are the series expansion analysis of Butera 
and Comi Q and by Campostrini et al. [10], and the field-theoretic analysis 
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of Guida and Zinn-Justin [11|. The latter work presents the result y^ = —0.81 
which is of some importance in scaling analyses of Monte Carlo data. 



2 The nearest-neighbor Ising model 

We consider the Ising model with spins s X; y :Z = ±1 on the sites (x, y, z) of 
the simple-cubic lattice, with nearest-neighbor couplings only: 

T~L/kT = —K nn ^ ' Sx,y,z{. s x+l,y,z s x,y+l,z + s x,y,z+l) (1) 
x,y,z 

The data used in this analysis of the nearest-neighbor Ising model are the 
result of Cluster Processor simulations (roughly one half processor-year), 
supplemented by some results for relatively small system sizes, obtained 
with workstations and PC's. These simulations used L x L x L systems 
with periodic boundaries and took place very close to the critical point, 
most of them at K nn = 0.2216545. The lengths of the Wolff simulations are 
shown in Table |]. This table presents the number of millions of samples 
#S taken per system size, and the number of Wolff steps #W before taking 
each new sample. One Wolff step means the formation and spin inversion of 
one cluster. These data are shown only for L > 40 in Table [l]. Also included 
in the present analysis are Monte Carlo data for some smaller system sizes, 
namely L = 4-16, 18, 20, 22, 24, 28 and 32. These data, as well as the 
details of the simulations, were already described in Ref. ||. 

Table 1: Length of Monte Carlo runs for 5 different Ising models, in millions 
of sampled configurations (#S). The number #W of Wolff clusters flipped 
before a new configuration was sampled is also indicated. Smaller system 
sizes L < 40 are also included in the present analysis (see text). 



K3 n /K nn 


0.0 


0.1 


0.2 


0.3 


0.4 


L 


#S 


#w 


#s #w 


#s 


#w 


#S 


#w 


#S #W 


40 


100 


10 


100 10 


100 


10 


100 


10 


100 10 


48 


50 


20 




50 


32 


50 


32 




64 


50 


32 




50 


32 


50 


32 




128 


18 


64 




20 


64 


20 


64 




128 


2.0 


100 














256 


2.1 


200 
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The 'data taking' during the simulations included the sampling of the 
magnetization m = L~ d J2 x ,y,x s x,y,x and the sum of nearest-neighbor prod- 
ucts (the energy of the nearest-neighbor model), and the accumulation of 
powers and cross-products of these. On this basis we can obtain a number 
of expectation values. Theoretical expressions for these quantities are also 
available. By matching these to the numerical data, one can verify the the- 
ory and solve for unknown parameters such as the location of the critical 
point. 

The numerical data include the dimensionless ratio 

Q L (K nn ) = ( 2 ) 

of moments of the magnetization m per spin in a system of size L. This ratio 



is a form of the Binder cumulant [12|, a quantity that has proven to be very 
useful in analyses of critical-point properties. We copy the renormalization 
prediction for the scaling behavior of QL{K nn ) from Ref. ||: 

Q L (K nn ) = Q + ai (K nn - K C )V* + a 2 (K nn - K c ) 2 L 2yt + 

a 3 (K nn - K c ) 3 L 3yt +••• + £ bjLV (3) 

3 

The sum on j is expected to be dominated by a contribution b\L yi where 
Hi = Vi ~ —0.83 [|l^, |3], 11] is the leading irrelevant exponent. In addition, 



finite-size scaling predicts terms [fj] 62 L V2 with 2/2 = 3 — 2yh ~ —1.963 and 
bsL y3 with 2/3 = y t — 2yh w —3.376. The term with exponent j/2 is due to the 
field-dependence of the analytic part of the free energy; a contribution with 
exponent 7/3 arises from nonlinear depencence of the magnetic scaling field on 
the physical magnetic field. Furthermore, one may expect higher order con- 
tributions from irrelevant fields; thus one may also include a term qib\L 2yi . 
Finite-size scaling also predicts further contributions to the sum in Eq. (^). 
Even if their exponents are smaller than y 2 or 2/3, this does not necessarily 
mean that they are less 'important' than 1/2 or 7/3, because the amplitudes 
may be very different. Which of these contributions are incorporated in 
the 'best fit' formula may be decided on the basis of such requirements as 
that the residual x 2 °f t ne least-squares analysis is acceptable, and that the 
fitted parameters become reasonably independent of the minimum system 
size cutoff L m i n of the data included in the fit. A considerable number of 
fits was tried out; in Table |2] we define 8 different types by specifying which 
of the parameters were included in the minimization procedure. In those 
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Table 2: Definition of several types of fits used in the analysis. The param- 
eters marked by '+' were actually solved in the minimization process; the 
remaining ones were kept at a constant value (zero except in the case of yi, 
see text). 



type 


Q 






62 


b 3 


y% 




1 


+ 


+ 


+ 










2 


+ 


+ 


+ 






+ 




3 


+ 


+ 


+ 


+ 








4 


+ 


+ 


+ 


+ 




+ 




5 


+ 


+ 


+ 


+ 


+ 






6 


+ 


+ 


+ 


+ 


+ 


+ 




7 


+ 


+ 


+ 


+ 




+ 


+ 


8 


+ 


+ 


+ 


+ 


+ 


+ 


+ 



cases where the irrelevant exponent yi was not fitted, it was set at the value 
Ui = —0.81 taken from Guida and Zinn-Justin |1 
The results of these fits are shown in Table 

2 it is necessary to discard an appreciable number of the smaller system 
sizes. It appears to be difficult to estimate the universal constant Q ac- 
curately; its variation between different fits exceeds some of the statistical 
error estimates. 

3 Including third-neighbor interactions 

We use the same reduced spin-one-half Hamiltonian as in Ref. || , namely 

x,y,z 

K3n\8x+l,y+l,z+l + S x +l,y-l,z-l + Sz-l.y+l^-l + (4) 

Whereas the analysis of Ref. j|] was restricted to only one nonzero value 
Ks a /K nn = 0.4, here we include similar Monte Carlo data for K% n /K na = 
0.1, 0.2 and 0.3 as well. For the smaller system sizes, the run lengths and 
further details of the simulations are mostly equal to those specified in Ref. 
H for K 3n /K mi = 0.4. For system sizes L > 32, the lengths of the Monte 
Carlo runs are included in Table |l|. The focus of these additional simulations 
was placed on the cases K^, n /K nn = 0.2 and 0.3, where the leading correction 
to scaling becomes small. 



]• 



For fits of type 1 and 
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Table 3: Results of least-squares fits to the dimensionless amplitude ratio 
Q, using the simulation data of the nearest-neighbor Ising model. Different 
types of fits are used; the type of fit is indicated in column 1. Data for system 
sizes below the specified minimum size (column 2) were ignored. Columns 
3-5 show the results for the critical point, the Binder cumulant and the 
irrelevant exponent. Where indicated by '(f)', y% was fixed at a constant 
value —0.81. The last two columns show the residuals and the number df 
of degrees of freedom of the fits respectively. 



type 


L > 


K c 


Q 


Vi 


x 2 


df 


1 


12 


0.22165434 


(5) 


0.62261 


(5) 


-0.81 


(f) 


54.4 


38 


1 


14 


0.22165438 


(5) 


0.62270 


(7) 


-0.81 


(f) 


45.6 


36 


1 


16 


0.22165442 


(5) 


0.62280 


(8) 


-0.81 


(f) 


38.8 


33 


2 


10 


0.22165462 


(6) 


0.6240 


(2) 


-0.95 


(2) 


43.1 


38 


2 


12 


0.22165456 


(7) 


0.6237 


(2) 


-0.91 


(2) 


35.7 


34 


2 


14 


0.22165457 


(8) 


0.6237 


(3) 


-0.92 


(3) 


35.2 


32 


3 


8 


0.22165451 


(6) 


0.62320 


(8) 


-0.81 


(f) 


50.3 


51 


4 


7 


0.22165459 


(7) 


0.6238 


(3) 


-0.91 


(4) 


55.0 


52 


4 


8 


0.22165453 


(8) 


0.6233 


(4) 


-0.83 


(6) 


50.2 


50 


5 


4 


0.22165455 


(5) 


0.62338 


(7) 


-0.81 


(f) 


56.1 


56 


5 


5 


0.22165457 


(6) 


0.62344 


(9) 


-0.81 


(f) 


54.7 


55 


5 


6 


0.22165456 


(6) 


0.62341 


(11) 


-0.81 


(f) 


54.5 


54 


6 


4 


0.22165456 


(8) 


0.6234 


(3) 


-0.82 


(5) 


56.0 


55 


6 


5 


0.22165450 


(9) 


0.6230 


(5) 


-0.74 


(7) 


53.7 


54 



Some fits of the types defined in Table § were applied to the Q data. 
The results were found to agree satisfactorily with universality of Q, so that, 
following Ref. ||, the combined Q data for the 5 models were analyzed by 
a single fit formula Eq. (||) in which the universal parameters appear only 
once, and the nonuniversal ones in 5-fold. This was done for several choices 
of combinations of parameters as defined in Table 0, in particular those 
types that include the irrelevant exponent y^ as a free parameter. Some of 
the results are shown in Table ||[ 

The simultaneous analysis of these 5 models is seen to allow a somewhat 
more precise determination of Q and y^. Again, the fit of type 2, involving 
only one correction to scaling, is seen to converge too slowly as a function 
of the minimum system size. The convergence is seen to improve when the 
number of parameters increases. However, the fits of types 7 and 8, which 
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Table 4: Results of least-squares fits for the dimensionless amplitude ratio 
Q, using the combined simulation data for 5 Ising models. Different types of 
fits are used (see text); the type is indicated in column 1. Data for system 
sizes below the specified minimum size (column 2) were ignored. The results 
for K c apply to the nearest-neighbor model. 



type 


L > 


K c 


Q 


Vi 


x 2 


df 


2 


14 


0.22165431 


(5) 


0.62234 


(8) 


-0 


780 


(9) 


148.2 


64 


2 


16 


0.22165437 


(5) 


0.62249 


(10) 


-0 


781 


(12) 


105.7 


53 


2 


18 


0.22165441 


(5) 


0.62267 


(12) 


-0 


787 


(16) 


79.0 


47 


4 


5 


0.22165444 


(5) 


0.62342 


(4) 


-0 


906 


(6) 


361.2 


164 


4 


6 


0.22165452 


(5) 


0.62352 


(5) 


-0 


887 


(8) 


221.4 


159 


4 


7 


0.22165454 


(5) 


0.62348 


(6) 


-0 


859 


(11) 


186.4 


153 


4 


8 


0.22165457 


(5) 


0.62348 


(7) 


-0 


839 


(15) 


175.6 


147 


6 


4 


0.22165462 


(5) 


0.62370 


(6) 


-0 


849 


(12) 


193.8 


164 


6 


5 


0.22165460 


(5) 


0.62358 


(8) 


-0 


824 


(17) 


184.8 


159 


6 


6 


0.22165458 


(5) 


0.62341 


(11) 


-0 


799 


(24) 


177.2 


154 


7 


5 


0.22165455 


(5) 


0.62367 


(5) 


-0 


815 


(7) 


213.3 


163 


7 


6 


0.22165459 


(5) 


0.62367 


(6) 


-0 


814 


(12) 


186.3 


158 


7 


7 


0.22165459 


(5) 


0.62357 


(7) 


-0 


797 


(19) 


176.0 


152 


8 


4 


0.22165463 


(5) 


0.62372 


(7) 


-0 


834 


(22) 


193.2 


163 


8 


5 


0.22165460 


(5) 


0.62358 


(9) 


-0 


826 


(37) 


184.5 


158 


8 


6 


0.22165459 


(6) 


0.62344 


(12) 


-0 


760 


(56) 


176.7 


153 



include a term q\b\l? Vl (not 5, only one because q\ is universal) in addition 
to those already given in Eq. ||, did not provide strong evidence for such a 
nonlinear contribution in the leading irrelevant scaling field. 

4 The magnetic and the temperature exponents 

A factor kT is included in the susceptibility, so that the high-temperature 
dependence vanishes and x = 1 according to Curie's law. It is thus expressed 
in magnetization fluctuations as 

X = L d (m 2 ) . (5) 

The expected finite-size scaling behavior is 

X = c + ci(iT nn - K c ) + ■ ■ ■ + L 2 ^- d [a + 
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Table 5: Results of least-squares fits to the finite-size data for the suscep- 
tibility of the nearest-neighbor Ising model. The type of fit is indicated in 
column 1. The minimum system sizes used in the fits are shown in column 
2. 



type 


L > 


Vh 


x 2 


df 


A 


5 


2.48222 


(17) 


52.0 


37 


A 


6 


2.48182 


(21) 


31.1 


36 


A 


7 


2.48164 


(25) 


27.6 


35 


A 


8 


2.48149 


(29) 


25.8 


34 


B 


4 


2.48151 


(25) 


30.0 


37 


B 


5 


2.48131 


(31) 


27.7 


36 


B 


6 


2.48139 


(38) 


27.4 


35 


C 


4 


2.48179 


(21) 


33.9 


37 


C 


5 


2.48148 


(27) 


26.9 


36 


C 


6 


2.48145 


(32) 


26.9 


35 



ai (K nn - K c )L y > + a 2 (K nn - K C )L 2 ^ + b x L^ + b 2 L* ■ ■ •] (6) 

Several types of fit were applied to the numerical susceptibility data for the 
3-d nearest-neighbor Ising model. In all fits, the critical point was fixed 
at K c = 0.22165459 and the irrelevant exponent at yi = —0.81. In fits of 
type A (see Table |||) we set b 2 = in Eq. @; in types 2 and 3, it was left 
free. The associated second correction exponent was fixed at yj = yt — 2yh 
(type B) or at yj = —2y^ (type C). The type B value follows when one 
includes a quadratic term in the physical magnetic field in the temperature- 
like scaling field; the type C value represents the field dependence of finite- 
size contribution to the critical free-energy density. 

As explained in Ref. || it is possible to obtain the temperature derivative 
of the Binder cumulant by sampling magnetization moments, the energy and 
their correlations. This quantity is expected to scale as 



9Q m 



u + Ul (K nn - K C )L^ + u 2 (K nn - K c fL 2 ^ + ■■■ + 



dK n „ 

uiLM + v 2 LW + •••]. (7) 

Several types of fit were applied to the Monte Carlo data obtained for this 
quantity in the case of the 3-d nearest-neighbor Ising model. The critical 
point was fixed at K c = 0.22165459 and the irrelevant exponent at y^ = 
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Table 6: Results of least-squares fits to the finite-size data for the tempera- 
ture derivative of the Binder cumulant of the nearest-neighbor Ising model. 
The type of fit is indicated in column 1. The minimum system sizes used in 
the fits are shown in column 2. 



type 


L > 


Vh 


x 1 


d f 


a 


5 


1 


5841 


(3) 


50.5 


43 


a 


6 


1 


5845 


(4) 


45.9 


42 


a 


7 


1 


5847 


(4) 


45.4 


41 


a 


8 


1 


5852 


(5) 


41.2 


40 


a 


9 


1 


5857 


(5) 


37.5 


39 


a 


10 


1 


5858 


(6) 


37.1 


38 


b 


4 


1 


5860 


(5) 


40.7 


43 


b 


5 


1 


5861 


(6) 


40.6 


42 


b 


6 


1 


5866 


(7) 


40.2 


41 


b 


7 


1 


5872 


(9) 


37.4 


40 


c 


4 


1 


5852 


(4) 


43.1 


43 


c 


5 


1 


5854 


(5) 


42.3 


42 


c 


6 


1 


5859 


(6) 


41.3 


41 


c 


7 


1 


5865 


(7) 


36.8 


40 


c 


8 


1 


5867 


(8) 


36.5 


39 


d 


4 


1 


5865 


(9) 


40.2 


42 



—0.81. In fits of type a (see Table ^) we set vi = in Eq. (^); in types b, c 
and d it was left free. The associated second correction exponent was fixed 
at yj = 3 — 2yh (type b) or at yj = y t — 2yh (type c). The type b value is due 
to the field-dependence of the analytic part of the free energy, and the type 
c value follows when one includes a quadratic term in the physical magnetic 
field in the temperature-like scaling field. Both terms were included in fits 
of type d. 

The fits of type a display a slow trend as a function of L so that ad- 
ditional corrections should be taken into account. The residuals do not 
clearly discriminate between type b or c corrections. The results for both 
types agree well with that for type d. 
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5 Discussion 



A major factor limiting the accuracy of Monte Carlo analyses of critical 
phenomena is the multitude of correction-to-scaling mechanisms that should 
be included in a proper theoretical description, i.e. in the fit formula that is 
used to model the finite-size data. Each correction that is neglected in the 
data analysis will lead to a bias in the results of the analysis. Unfortunately, 
the statistical accuracy of the estimated parameters tends to decrease rapidly 
when more free parameters are introduced. It is thus necessary to restrict 
the number of free parameters. One may estimate the effect of the neglected 
corrections by imposing a minimim system size L m i n , ignoring data for L < 
L min , and observing the way in which the residual x 2 depends on L min . 
Fortunately, x 2 tends to decrease rapidly when L m \ n increases, especially 
when the number of parameters is not too small. 

Very few free parameters appear in the fits of types 1 and 2 in Tables ^ 
and ^| only one finite-size correction is taken into acount. This is reflected 
by the large values of L m \ n needed for acceptable fits, and by the large 
differences of e.g. the irrelevant exponent y^ with respect to the expected 
value ~ —0.82. Similar effects are seen in the recent analysis of Ballesteros 
et al. who also use one correction. The effects on our present analysis are 
more serious because the data have a somewhat higher statistical accuracy. 

Thus, we found it necessary to include a second correction in the analysis 
of the Binder cumulant, the term with exponent 2/2 = 3 — 2yh in Eq. ([|). 
This leads to a large reduction of the residual x 2 \ the fits (of types 3 and 
4) are now seen to become acceptable already at L m \ n = 8 and then yield 
reasonable values for j/j. 

It is remarkable that the term with exponent y% was not present in the 
analysis given by Hasenbusch et al. Q, which nevertheless included small 
system sizes. This term is likely to be responsible for the 'slow convergence 
of the NNN flow' observed by Hasenbusch et al. in their analysis of an 
Ising model with first and third neighbor interactions. No such ill effect was 
however observed when the same analysis was applied to a spin-one model 
||. The different behavior of these two models is in agreement with results 
already given in Ref. ||] which pertain to models that do not differ much 
from those of Ref. [|| . These results include the finite-size amplitudes of the 
L V2 term (62 in Tables 5 and 6 of Ref. || ) ; this amplitude is appreciable for 
the Ising model with first and third neighbor interactions, and approximately 
zero for the spin-one model. 

The present results in Tables ^ and |] illustrate that a further strong 
reduction of x 2 occurs when a third correction term is included. The fits of 
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Table 7: Final results of the present work. The argument of the quantities 
given in the first column is the coupling ratio K^ n /K nii . The values of the 
critical points K c apply to the nearest-neighbor coupling K nn . The ampli- 
tude b\ of the leading finite-size correction in Q was determined assuming a 
fixed irrelevant exponent yi = —0.81. Two-sigma error estimates are given 
in the rightmost column. 



quantity 


value 


error 


K c {0.0) 


0.22165459 


(10) 


K c (Q.l) 


0.18562466 


(52) 


K c (0.2) 


0.16073235 


(12) 


Kjo.3) 


0.14230187 


(12) 


K c (0A) 


0.12800393 


(40) 


6i(0.0) 


0.0896 


(25) 


6i(0.1) 


0.0494 


(25) 


6i(0.2) 


0.0142 


(25) 


6i(0.3) 


-0.0162 


(25) 


6i(0.4) 


-0.0465 


(25) 


Q 


0.62358 


(15) 


Vh 


2.4814 


(5) 


Vt 


1.5865 


(14) 


Vi 


-0.82 


(3) 



types 5, 6 and 8 are able to describe the finite-size dependence of Q for even 
smaller values of L. It is not completely clear that the term with exponent 
U3 = Vt — tyh is actually the most important that has to be accounted for. 
There may also be significant corrections due to a second irrelevant field 
with an unknown exponent, and due to the expected nonlinear dependence 
of the magnetic scaling field on the physical field. The latter mechanism 
produces a finite-size correction in Q with exponent = —2y^ which is 
even smaller than y%. Fortunately the important fitted parameters do not 
depend strongly on the choice between y% and y±. This observation gives 
us some more confidence that the statistical error estimates are reasonable 
indicators of the actual uncertainties. We still prefer to give conservative 
two-sigma error bounds, in order to allow for the arbitrariness in the form of 
the fit formulas. Our preferred values and error estimates are summarized 
m Table 0. 

The result for the critical point of the nearest-neighbor model agrees well 
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with that of Ballesteros et al. |7j and reasonably with that of Hasenbusch et 
al. H . Also the universal quantities given in Table [7| are mostly consistent 
with Refs. M and M but we note a rather large difference with Q of Ref. 
H]. A possible explanation of this difference is the presence of the term 
62 L y2 in Eq. (||) which was not taken into account in Ref. B. Although the 
amplitude 62 is small ||, we see no reason why it should precisely be zero. 

Although our present work involves a body of simulations that is larger 
than those of Refs. J7|] and ||, this is not reflected in a smaller error bar 
for yh (see Table 0) than those given in Refs. § and g. The mam reason 
is obviously that our quoted errors are two standard deviations. The errors 
given in Refs. @ and || should be multiplied by two before comparing with 
our work. 

Furthermore we note that our results for the critical exponents agree 
well with recent series-expansion results of Campostrini et al. [10]. Finally, 
our results for yi and yt are in a good agreement with the field-theory results 



obtained recently by Guida and Zinn- Justin [11]. Although the agreement 
is not so good in the case of y^ where the difference is slightly larger than 
the combined error bars, the overall agreement between the results of such 
different approaches is quite satisfactory and confirms the underlying as- 
sumptions of scaling and universality. 
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